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Abstract. We study the primordial magnetic field generated by the simple model f^FF 
in Starobinsky, i?^-inflationary, model. The scale invariant PMF is achieved at relatively 
high power index of the coupling function, \a\ ^ 7.44. This model does not suffer from the 
backreaction problem as long as, the rate of inflationary expansion, i7, is in the order of 
or less than the upper bound reported by Planck (< 3.6 x 10“^Mpi) in both de Sitter and 
power law expansion, which show similar results. We calculate the lower limit of the reheating 
parameter, R^ad > 6.888 in i?^-inflation. Based on the upper limit obtained from CMB, we 
find that the upper limits of magnetic field and reheating energy density as, (p5end)cMB ^ 
1.184 X and (preh)cMB ^ 8.480 x 10“^^Mpp However, the limits derived from 

the inflationary model are (pSend)R2-inflation < 4-6788 x and (Preh)R2_inflation < 

3.344 X 10 ^^Mpp All of foregoing results are well more than the lower limit derived from 
WMAP7 for both large and small field inflation. By using the Planck inflationary constraints, 
2015 in the context of i?^-inflation, the upper limit of reheating temperature and energy 
density for all possible values of are respectively constrained as, T^eh < 4.32 x lO^^GeV 
and preh < 3.259 x 10“^^Mpj at 77/g 0.9674. This value of spectral index is well consistent 

with Planck, 2015 results. Adopting Treh, enables us to constrain the reheating e-folds 
number, A^reh fhe range 1 < A^reh < 8.3, for —1/3 < iCreh < 1- By using the scale 
invariant PMF generated by pFF^ we find that the upper limit of present magnetic field, 
Ho < 8.058 X 10“^G. It is in the same order of PMF, reported by Planck, 2015. 
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1 Introduction 

Inflationary cosmology solves the fundamental problems of the Big Bang model [1]. It pro¬ 
vides an explanation for the large scale structure of the Universe, which is linked to the 
quantum fluctuation of the field of inflation, 0. The implications of inflation have been 
consistent with the observational aspects of the Universe in the large scale. 

On the other hand, large scale magnetic fields are being observed in all kinds of galaxies 
and cluster of galaxies at wide range of redshifts. Some evidences for the presence of the 
magnetic fields in voids were reported [2]-[8]. Such detection strongly enhances the primor¬ 
dial origin of the magnetic fields (PMF). The latest constraints on the value of PMF were 
presented recently by Planck 2015 [9]. Hence, the seed of magnetic field cannot be generated 
by the galactic dynamo in the absence of charged particles [10]-[12]. The generation of PMF 
is still an open question in cosmology, see the reviews of this subject [13]-[20]. However, 
one of the most interesting cosmological models is the simple (/^FF)model [21]-[25]. It has 
gained more interest because it is a stable model under perturbation. (See [26] and references 
therein.) Also, it can lead to a scale invariant spectrum of PMF [23]-[25]. A scale invariance 
property explains why PMF is detected nearly in all scales of the universe. 

The main problems with this model are: the backreaction problem, where the scale of 
the energy of the associated electric field with PMF can exceed the scale of inflation itself 
[24]-[28], and the strong coupling between electromagnetic fields and charged matter at the 
beginning of inflation [27]-[28]. There are some computational and phenomenological ways, 
which were proposed to avoid these problems [22]-[28]. For example, as the amplitude of 
magnetic field spectrum reach around the scale invariant value (^ nG) during inflation, the 
total electromagnetic energy density is much less than the energy of inflation itself [22]. 
However, the avenues, proposed to overcome these two problems, need more investigations in 
order to form a robust model, which is self-consistent in both classical and quantum regimes. 

The first problem is classical and may spoil the inflation. However, the second problem 
is of quantum nature and may lead to a huge coupling between the electromagnetism and the 
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charged particles [27]. Thus, if we define this model as a classical consistent one, it cannot 
be ruled out only because of a quantum effect like strong coupling. Nevertheless, proving 
the consistency between the quantum limit and the classical limit of the model will make it 
more trustable. For this reason, solving the classical problem (backreaction) associated with 
this model makes it self-consistent model. That is what we try to achieve as a result of this 
paper. 

The Lagrangian of a scalar (inflaton) field (j) coupled to the gauged electromagnetic 
vector field can be written [23]-[25] as, 

^ = - Q ^ (0) + (0, t) (1.1) 

where, — djyA^ — d^Ajy is the electromagnetic field tensor and g is the determinant of 
the spacetime metric g^y. The first term in the Lagrangian is the standard kinetic part of 
the scalar field, and the second term,!/ (0), is the potential, which specifies the model of 
inflation. 

The Lagrangian of a pure electromagnetic field is of the form, — In Eq.(l.l), we 

couple it to the scalar field through the unspecified function /((/), t). The main reason behind 
this coupling is to break the conformal symmetry of electromagnetism and hence prevent the 
dilution of the seed of magnetic field as it is generated in the inflation era. On other word, 
electromagnetic fields would not feel inflation. At the end of inflation, the coupling function, 
/((/), t) ^ 1, to retrieve the conformal electromagnetism. The last condition is important to 
decide the form of coupling function and to study PMF in post-inflationary phases. 

The standard model of inflation based on a single scalar field, such as quadratic,!/ (0) ^ 

0^, quartic, V {(j)) ^ 0^[29], i?^-inflation, V (0) ^ ^1 — exp [—^y2/3(|)/MpL]^ [30, 41] and 

the exponential potential, V (0) ^ exp[— v^2ei (0 — (/)o)][31]. The last one is used in [23]-[24] 
to find the magnetic and electric spectrum in the f^FF model. These models became more 
interesting after WMAP [32]-[33] and Planck [34]-[35]. 

Very recently, the joint analysis of BICEP2/Keck Array and Planck (BKP) data was 
released on Feb 2015. The joint data of three probes eliminate the effect of dust contamination 
and show that the upper limit of tensor to scalar ration, ro.os < 0.12 at 95% CL, and the 
gravitational lensing B-modes (not the primordial tensor) are detected in 7a [36]. Similarly, 
the inflationary models are constrained based on the new analysis of data. The scalar spectral 
index was constrained by Planck, 2015 to be rig = 0.9682 ±0.0062 [37]. As a result, the more 
standard inflationary models, like i?^-inflation, which result low value of r, are the most 
favored one. However, the chaotic inflationary models like large field inflation (LFI) and 
natural inflation (NI) are disfavored [37]. These results ruled out the first results of BICEP2, 
2014 [38]. 

In this paper, the simple model, f^FF^ of PMF will be investigated in the context of 
i?^-inflation, in the same way we did for NI and LFI, [39]-[40]. Further, we will constrain the 
reheating parameters under the same inflationary model by using the reported upper limit 
of PMF by Planck, 2015. The present PMF will be constrained based on the scale invariant 
magnetic field generated during inflationary era. 

Throughout this paper, we adopt the natural units, [c^h^kB — 1], the signature 
( — 1,1,1,!), and flat universe, where we use the reduced Planck mass, Mpi = (87rG)~^^^. 
It will be taken, Mp\ = 1 in the computation parts. Hence, the potential of i?^-inflation 
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(Starobinsky model) can be written in Einstein frame [41]-[42] as, 

^(0) =M4(l-exp[-y^</./Mpi])^ (1.2) 

where, M is the amplitude of the potential and it can be determined by the amplitude of 
CMB anisotropies. 

The order of this paper will be as follows, in section. 2, the slow roll inflation formulation 
is presented for both simple de Sitter model of expansion and the more general power law 
expansion in the context of i?^-inflation. In section. 3, the PMF and associated electric fields 
are computed. The constraining of reheating parameters in the same model of inflation and 
by using the reported results of Planck, 2015 are discussed in section. 4. In section. 5, the 
present PMF is constrained based on the computed magnetogensis. In section. 6, a summary 
of the results is presented. 


2 Slow roll analysis of i?^-inflation 


During inflation, we will assume the electromagnetic field to be negligible compared to the 
scalar field, (j) [24]. Hence, the equation of motion derived from 1.1 for the scalar field can be 
written as. 


(j) + 3H(j) + — 0 


( 2 . 1 ) 


where, H (t) = a{t) /a{t)^ is the Hubble parameter as a function of cosmic time, t, and 
a (t) is the cosmological scale factor. The over dot indicates differentiation with respect to 
cosmic time, and The Friedman equation can be obtained from the Einstein field 

equations by assuming a spatially-flat Friedmann-Robertson-Walker FRW universe, which 
yields. 
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where, is the curvature of the universe and. A, is the cosmological constant. The last two 
terms of 2.2 can be neglected in the inflation era. Also, under slow roll approximation, one 
can neglect the second derivative in (2.1), which leads to the attractor condition. 
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They also, can be written in terms of the Hnbble parameter, 
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Finally, the relation between the two formalisms [43]-[44] can be written as 


^iv = 


/ 3 - 62 H 
\3 - eiH 


2 


All of the above parameters are assumed to be very small during the slow roll inflation, (eiy, 
and 62 if) <C 1. Further, inflation ends when, the values of (eiy, ein) 1. In the first 
order of approximation, one can neglect ein and e 2 ifcomparing with 3, obtaining ew — cih- 
Therefore, using (2.7) and the relation between the cosmic time, t, and the conformal time, 77 , 
dt — a ( 77 ) df]^ one can write the relation between conformal time and slow roll parameter, 
eiH, as [45], 


Assuming, cih ^ const^ and then integrating (2.9) yields the power law expansion of the 
universe during inflation, 

a{ri) = , ( 2 . 10 ) 

where, I is integration constant. 

On the other hand, the universe expands exponentially during inflation at a very high 
but almost constant rate, i7, in the simplest form of inflationary expansion (de Sitter), 

= - ~ const, ( 2 . 11 ) 

a 


a{t) — a {ti) exp [Hit] , ( 2 . 12 ) 

where, ti, is the start time of inflation. In a conformal time, Eq.(2.12) can be written as, 

«(??) =- 77 -Vv’ 

H T] + Cl 

where, ci, is the integration constant. Plugging ([15]) into the relation between cosmic and 
conformal time and integrating, implies that 77 ^ (— (X),0“) as t ^ (0,(X)). Thus, if 77 ^ 0“, 
then a ( 77 ) ^ (X) and ci ^ 0. Therefore, we can write a ( 77 ) = — 1/ (H 77 ). 

Also, the relation between slow roll parameters and the scalar power spectrum amplitude,A^, 
the tensor power spectrum amplitude,At, the scalar spectral index, 77 . 5 , and tensor-to-scalar 
ratio, r, can be written as follows [43]-[46], 
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Figure 1. The r — Ug relation in i?^-infiation for N = 50,60,64, and 68. All values of r^v at 
Us — 0.968, are well below upper limit of r(< 0.12) by BKP, 2015. 


One can find the relation between r and ris which depends on the number of e-folds of 
inflation, N. The first order of approximation for N can be written as, 


A^ = ln 


a{v) 

aiVi) 
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0end 1 3 

2 / — 7 


2Mpi 0 


2 (j) 


3 Afpi 


(2.18) 


where, 0, is the initial field and (^end is the field at the end of inflation. Solving for 0 from 
(2.18), 




(2.19) 


where, A^end is the number of e-folds at the end of inflation. For simplicity, we will denote 
(A'end — N) = N. Hence, plugging (2.19) into (2.16)-(2.17), yields. 
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( 2 . 20 ) 


The relation (2.20) is drawn for some interesting values of N in Fig.l. The values of r, are 
well below the limit of BKP, 2015 (r < 0.12). 

3 The electromagnetic spectra generated in i?^-inflation 

One can write the electromagnetic equation, for d^from the Lagrangian, (1.1), as. 


da 




= 0 , 


(3.1) 


where, g is the determinant of the flat FRW space-time metric defined as, 

ds^ — gfjiydx^ dx^ = —dt^ + c? (t) dx^ = dx^ = —a^ {g) {dg^ -h . (3.2) 
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Hence, (t), y^—g' = {rf)^ and — S^. Adopting the Coulomb (radiation) 

gauge, diA^ (t, x) = 0, and the charge-free condition Aq (t, x) = 0, Eq.(3.1) can be written as 


Ai{t,x) + (h +“^1 Ai {t, x) - djd^Ai {t, x) = 0. 


(3.3) 


In conformal time, the above equation can be written as 

A"i {ri, x) + 2 y A'j ( 77 , x) - ( 77 ) djd^Ai {rj, x) = 0 . 

Define the function, ( 77 , x) = f ( 77 ) ( 77 , x), and plug it into (3.4), 

A"j + ~ dj&’Ai (rj, x) = 0 . 


(3.4) 


(3.5) 


On the other hand, the quantization of A^ can be written in terms of creation and 
annihilation operators, and hx (A:), as. 


r ^ 

Aj (? 7 , x)= -—^ ^ eix {k) [ 6 a (k) A ( 77 , k) + ^^^a {k) A* ( 77 , k) (3.6) 

J (27r)"/“ 

where, Cix is the transverse polarization vector, and A: = ^, is the commoving wave number. 
Substituting (3.6) into (3.5), 


A" ( 77 , fc) + ( - y ) ^ ( 77 , k) = 0 , 


(3.7) 


where we have used, = g^^dk = Similarly, the stress energy tensor, can be 

written in terms of the action, as, 

= = {d^4>) {d,4>) - f {4>,t)FfS5Fa^. (3.8) 

For the magnetic field, we ignore the kinetic part (first term) because it does not contribute 
to the electromagnetic field. Then magnetic part of 

T^o = (</>, t) (djAi - diAjXdiAk - dkAi) = (3.9) 

The energy density of the magnetic field, pb is found by taking the average of the 
component of the stress energy tensor. Finally, the spectrum of the magnetic field is given 
[24] by 


dpB _ J_ f 


fc|^(7?,fc)| . 


dink 277^ \ ® 7 
Similarly, the spectrum of the electric field is given by 


dpE f k^ 
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Therefore, to calculate the electromagnetic spectra, one has to find ^( 77 ,/c), which can be 
achieved by solving (3.7). But first, we need to specify the coupling function, / (0). Assuming 
that the relation between the coupling function and scale factor is of the power law form [ 8 ], 
/ ( 77 ) oc where, a is free index to be determined later from the shape of the spectrum 
of PMF. Then, by combining (2.1) and (2.2) in the slow roll limit. 


/ {(/)) oc exp 


a 


0 V 

3 Mpi^ ^ 


(3.12) 


Substituting (2.19) into ( 1 . 2 ) and (3.12) yields the coupling function as a function of A^, 


f{N) =D 



(3.13) 


where, D is a, coupling constant. Substituting (3.13) into (3.7) gives, 

(r?, k) + {e-Y (r?)) A (rj, k) = 0 , (3.14) 

here the function Y ( 77 ) = Hence, we need to write the derivative in terms of e-folds 
number, N. In the next two sections, the electromagnetic spectra are calculated for both de 
Sitter and power law expansion. 


3.1 Inflationary electromagnetic spectra in de Sitter expansion 

Shortly after the onset of inflation the value of H becomes very high and is approximately 
constant, but later on, it decreases as the value of the field changes. For the zeroth approxima¬ 
tion and after the first few e-foldings, we can consider H as a constant ratio. That is basically 
the de Sitter expansion, which is exactly exponential expansion as described by (2.12). In 
fact, de Sitter model does not have graceful exit from inflation [46]. But it can be assumed 
as a valid approximation on part of the inflation. Assuming, a{r]i) = const, and substitut¬ 
ing of (2.13) into (2.18) and differentiating both sides, yields that, drjr]f{v) — V~‘^9NNf{N). 
Therefore, the function Y ( 77 ) can be written as. 


Yds{v) = 


a 


(3 - ANfa - 36 

144/VV 


(3.15) 


Also, as the variation in N is very small comparison with the variation in 77 , we will assume 
that N is quasi-constant and we will not write it as a function of conformal time explicitly 
in Eq.(36). Therefore, substituting of (3.15) into (3.14), yields. 


A iv, k) = (fcr?)'/2 [Cl ik) {krj) + C 2 (k) (fcr?)], (3.16) 


where, {kp)^ is the Bessel function of the first kind, and the argument of the function. 


Xds{(x,N) = 


12N 


(3.17) 


The relevant PMF is obtained in the long wavelength regime. A: <C 1 (outside Hubble 
radius). In this limit, Eq.(3.16) can be written as [24], 


A<1 {r], k) = (k) \Di (x) (/ct/)^+ 2 + L>2 (x) (kv) ^ 


2-x 


(3.18) 
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Figure 2. The a — N relation, at x(A/, a) = 5/2, the case at which we can generate the scale 
invariant PMF. 


The constants, Di (x) and D 2 (x)? can be fixed by using the normalization of A{r]^k) and 
the other limit, 


A >1 (ri, k) ^ e 


. Thus, they can be written as 

g-i7r(x+|)/2 




2 x r(x + l)cos(7r(x+ 5 )) 


By substituting (3.18) into (3.10)-(3.11), we get the spectra of the magnetic and electric 
fields. If the PMF is scale invariant, then the magnetic spectrum should be constant, which 
could be achieved at x = 5/2. This value can only be obtained around, a ^ {—7.44, 7.44}, 
see Fig. 2. 

Although these values of a are too high, which may exceed dynamo limit [24], we adopt 
the positive value, a ^ 7.44. In the last section, we will discuss this point. So, let us assume 
it is valid at this point and use it to investigate the electromagnetic spectra generated by 
f^FF for long wavelength approximation {k <C 1). From (3.10)-(3.11), one can plot the 
electromagnetic spectra for different values of H and cr, see Fig. 3, for H = 10“^Mpi, which 
is around the Planck, 2015 upper limit of pivot Hubble parameter, < 3.6 X 10“^Mpi) 
[37]. The pivot moment is the time when the commoving scale (fc* = 2 :^ 0.05Mpc“^) 

exits the Hubble radius. In this case, the scale invariant PMF can be generated and the 
backreaction problem may be avoided at cr ^ 7.4359. However, for a — 2^ the backreaction 
problem can be avoided easily, but the scale invariant PMF is not generated, see Fig. 3. 
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Figure 3. The EM spectra for a = [2,7.4359], A/" = 64, 77 = —10“^ and the expansion rate , 
H = 10“^Mpi. For a = 2, the scale invariant PMF cannot be achieved although the backreaction 
problem is easily to be avoided. However, for x(A^, a) = 5/2 {a = 7.4359), the scale invariant PMF 
can be generated and the magnetic field energy is less than the electric field which can avoid the 
backreaction problem. 


On the other hand, in the large expansion rate case, H = 0.2Mpi, the scale invariant 
PMF can be generated at a ~ 7.4359 but the problem of backreaction may not be avoided 
in the low, A; <C 1, see Fig. 4. In fact, this value of H is way above the constraint upper 
limit [37], but we use it to decide the scale at which the electric field exceeds magnetic field 
{k ^ 1.5 X 10“^). Hence, we can use this value of the wavenumber as an upper limit in 
this case. Again, for a = 2, the backreaction problem can be avoided easily, but the scale 
invariant PMF is not maintained throughout the inflationary era. 

The value of the index, cr, at which the scale invariant PMF can be generated, varies as 
the value of N changes. There is a slit different between this relation in de Sitter and power 
law expansion, see Fig. 5. These values are around, a ^ 7.44, for an interesting values of 
e-folds, (50 < N < 70). The validity of a ^ 7.44 will be discussed in the section. 6, but from 
now onward, we adopt this value at = 64 to study the scale invariant PMF. In general, 
the whole results of de Sitter are very close to the power law results but both of them will 
be presented for a completeness. 

First, let us draw the relation between electromagnetic spectra, Hubble parameter, i7, 
and e-folding number, A^, see Fig. 6. It is clear, that, for the values of H < 0.2Mpi, the 
backreaction problem may be avoided for all interesting values of N. Therefore, one can 
consider H ^ 0.2Mpi, as an upper limit of expansion rate below which, the backreaction 
problem can be avoided. The current investigation will be bounded in 10“^Mpi < H < 
0.2Mp\. 

Similarly, drawing the relation between the electromagnetic spectra and e-folding num¬ 
ber around the upper limit, H ^ 0.2Mpi, shows that below the value, N ^ 68, the backreac¬ 
tion problem can be avoided, see Fig. 7. Therefore, this value of A^ ^ 68, can be considered 


- 9 - 



^ ~ 7.4359 



1.5x10- 



10 “^ 


10 


-4 


k 


Figure 4. The EM spectra for a = [2,7.4359], = 64, 77 = —10“^ and the expansion rate, 

H = 0.2Mpi. For a = 2, the scale invariant PMF cannot be achieved although the backreaction 
problem is easily to be avoided. However, for x(A^, a) = 5/2 {a ^ 7.4359), the scale invariant PMF 
can be generated but the magnetic field energy is less than the electric field for low wavenumber. The 
PMF exceeds the electric fields at, /c ^ 1.5 x 10“^. 


as an upper limit below which, a scale invariant PMF does suffer from the backreaction 
problem. 

3.2 Inflationary electromagnetic spectra in power law expansion 

For more optimal slow roll analysis that has a smooth exit from inflation, one can use the 
power law expansion at which the Hubble parameter is not constant but a function of field, 
If the field falls below a certain value, it starts to oscillate and then converts into 
particles in the reheating era, right after inflation. The scale factor follows the power law, 
(2.10). In the slow roll limit, we can approximate . Hence, substituting (2.4) into (2.10) 
yields, 


_ ^ _4_ 

a{r]) = lr] (3.20) 

where, I is the integration constant. For relatively high (N > 50), the absolute value of 
the exponent in (3.20) is < 1.0003, so we can approximate I H~^. Therefore, one cannot 
expect significant differences between the two cases. Nevertheless, we will investigate this 
case to have more precise results. 

Again, as we did in the last section, we will assume that N is quasi-constant and we will 
not write it as a function of r] explicitly. Hence in the power law expansion, we substitute 
(2.10) into (2.18) and differentiating both sides, yields that, 


Ypl (r?) 


(3 - 4Ar)2(63 - 4Ar (63 + AN (-9 + AN))fa (-36 + (3 - ANfo^ 
144iV2(21 + SN (-3 + 2iV)) V 


(3.21) 
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Figure 5. The relation between a and N at which x(A/', a) = 5/2, for both de Sitter and power law 
expansion. For the interesting range of A^(50 <N< 70), the value of a is around, a ^ 7.44.. 
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Figure 6. The relation between the EM spectra, Hubble parameter i7, and e-folds number, N for 
de Sitter expansion. Around H ^ 0.2Mpi, electric field can exceed the magnetic field in all interesting 
range of ^"(50 <N< 70). Hence, one can consider H ^ 0.2Mpi, as an upper limit under which the 
backreaction problem can be avoided. 
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Figure 7. The relation between the EM spectra and e-folds number, N for both de Sitter and power 
law expansion around the value H ^ 0.2Mpi. The electric field can exceed the magnetic field around 
N 6S. Hence, one can consider it as an upper limit of e-folds number before the end of inflation 
under which the backreaction problem can be avoided. 
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Figure 8. The EM spectra for a = 7.4345, = 64, 77 = —10 ^ and the expansion rate, H ^ I ^ = 

0.2Mpi. The PME exceeds the electric fields at /c ^ 1.5 x 10“^. 


Plugging (3.21) into (3.14) yields the same solution as (3.16). In this case, Xpl{(^^N) = 5/2, 
at (a {—7.4339,7.4345}, see Fig. 2 for all interesting values of N. Therefore, the scale 
invariant PME can be obtained in the long wavelength regime, A: <C 1 (outside Hubble 
radius). Also, the electromagnetic spectra can be obtained in the same way as done in last 
section. The results is very close to the those obtained by de Sitter case, see Fig. 8-9. 

It is clear from Fig. 9 that, at = 3.6 x lO^^Mpp In fact that 

value of the Hubble parameter is the upper value during inflation or at the time of pivot 
scale, at which the space time exits the Hubble radius [37]. Therefore, the problem of the 
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Figure 9. The EM spectra for ol = 7.4345, N — 64, 77 = —10“^ and the expansion rate at the onset 
of inflation, = 3.6 x 10“^Mpi. The relation, stays valid, throughout inflation 

at low values of, /c ^ 1 . 


—|,dS&PL, jt= lO-'O W =64,()'^ 7.435.H^l-'=3.6x lO-^, ;c« 1 
dink dink 

m-lk, _ 


3.842x1 0 -'’m4 


Power Law 
-de Sitter 



Figure 10. The EM spectra for a = 7.4345, = 64, /c = 10“^^ and the upper limit of Hubble 

parameter during inflation, H ^ = 3.6 x 10“^. The relation, ^ stays valid, throughout 

the period, r] —1 . 


backreaction caused by the divergence of the electric field at A: <C 1, is avoided. In order to 
make sure, the above relation stays valid throughout the inflationary era, we can plot the 
EM spectra versus the conformal time, at 77 <C — 1, see Fig. 10. 

So far, we make sure that the electric field has negligible contribution during the infla¬ 
tion. However, to test the contribution of the magnetic field in the inflationary energy, one 
has to consider the constraints of post-inflation eras. That point will be discussed in the next 
section. 

4 Constraining reheating parameters by PMF 

It is widely believed, that there was a phase of pre-heating or reheating at the end of inflation 
and before the radiation dominated era [47]-[49]. In this phase, as the temperature of infla- 










tion falls to certain value, T^eh, inflaton starts oscillates around some value and decays into 
standard matter, that populates the Universe later. As the temperature continues to fall, 
the Big Bang Nucleosynthesis (BBN) starts taking place at Tbbn ^ IMeV. These produced 
particles are perturbatively decayed into radiation in the radiation era. 

Initially, this process was thought it only occurs in a complicated inflationary model that 
has more than one field [50]-[52]. However, later on it was shown that it could be occur even 
in a single field model but at sub-Hubble scale of perturbations [53]-[55]. In order to constrain 
this phase, one has to define the reheating parameters based on both the inflationary model 
and the observations, like the Cosmic Microwaves Background (CMB) [56], the Large Scale 
Structure (LSS) [57], the BBN [58], and magnetogenesis [59]. 

There are two main difficulties in investigating this era. First, no direct cosmological 
observations can constrains the reheating parameters. Second, the physics of this period 
is highly uncertain. Therefore, indirect constraints of the reheating parameters are usually 
calculated from other cosmological observations. Also, several models have been proposed 
to explain this era. In this section, the effect of the scale invariant PMF on the parameters 
of reheating is investigated. That is basically similar to Ref [59] but by considering R^- 
inflationary model. We adopt the new upper limits of present PMF which was constrained 
by Planck, 2015 [9] and the instantaneous transition to the reheating and the epochs come 
after. For this reason, we start by discussing the reheating parameters. Next, Next, we 
discuss how to constrain them by i?^-inflation and the present upper limit of PMF. 


4.1 Reheating analysis 

The reheating era can be specified mainly by three parameters, the reheating parameter, Rrad? 
the reheating temperature, T^eh, and the equation of state parameter, rcreh [43],[60]-[61]. The 
first one is defined as, 

/ \ 1/4 

l^end I Pend \ 


^rad — 


(4.1) 


\ / ’ 

l^reh V Preh / 

where, p is the energy density, and “end” and “reh” stand respectively to the end of inflation 
and the end of reheating era. From the conservation of energy during the reheating era, 
P = P 0 + P 7 and P — P(i) + P 7 / 3 , where is the energy density of the scalar field of inflation 
and p^ is the energy density of the radiation. Assuming these are the main constituent of 
the reheating era, one can write. 


Preh — Pend GXp 




Nreh _|_ w{a)] 


dCL / — Pend6^P[ 3AA^(1 + lypgh)]) 


'A^end “ 

where, lUreh is mean value of the equation of state parameter {P/p = w{a)) defined by, 

1 


(4.2) 




/ 

Jn, 


-da. 


A^reh A^end J ^ 

Taking the logarithm of (4.1) and making use of (4.2)yields, 

1 (Apeh ~ ^ \ 1 ~ -i f Preh A 

In i?rad = -^- (-1 + 3^reh) = -z-r In ' ' 


(4.3) 


12(1 -|- Wj;e\i) VPend 




In terms of the pivot quantities [30] , 
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1/4 

Pi 


7 . 


3I4nd 3 - eiy* ^ 


3 ^IVend 
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On the other hand, the relation between the reheating temperature and energy density 
can be written as, 


TT 


Preh — 2Q^reh?leh5 


(4.6) 


where, is the number of relativistic degree of freedom at the end of reheating. One also 
can relate the reheating temperature to today temperature of CMB, Tq, [62] as 


Treh = To 




V ^reh / V ^ f ^reh 


43 ^ 


^eq / \ 1 l^reh / 


1/3 


(4.7) 


where, aeq/areh = during radiation era, and “eq” stands for the period of equality 

between radiation and matter dominant phases. 

In terms of the pivot scale, k = dkHk^ at which the commoving scale crosses the Hubble 
radius (horizon size) during the inflation, one can write down the relation between the scale 
factors and Hubble parameters [43] for different epochs of the Universe as 


aoi^o cloHq 

From (4.8), one can write the ration. 


^ _ Cik ^end ^reh ^eq Hk 

^end ^reh ^eq ^0 Hq 


ao 


^eq 


g-iVfc g-iV,eh g-iVeq _ 

k 


(4.8) 


(4.9) 


In (4.7) and (4.9), we assume a form of exponential expansion takes place during reheating 
and radiation epoch too. Substituting of (4.9) into (4.7) gives. 


Tleh = 


(w 


J3 V/=>„ 


-^fcg-^reh 


(4.10) 


If the equation of state is assumed to be constant during reheating then, p oc a ^(1+^). 
Also, the relation between the energy density and the potential [63] can be written as. 


P= 1 + 


1 


3/eiy - 1 


U(0). 


(4.11) 


Making use of these relations and substituting of (4.10) into (4.6) and taking natural loga¬ 
rithm yields the equations of reheating e-folds numbers and temperature. 
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\ Vnd 

1 -| 

3(l+^reh) 
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3(1+Wreh) 

3^reh-l 


, (4.13) 


where, V^nd is the potential at the end of inflation and Ag is the scalar power spectrum 
magnitude, obtained by (2.14). Hence, from (1.2) and the condition, Cnv ~ 1? af fhe end of 
inflation, one can decide V^nd- In fnct, the field (/)end hns three different values, based on which 
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slow roll parameter is equal to unity. They are respectively, (/)end — 0.94Mpi, 1.83Mpi, l.SlMpi 
for €iy, 62 v stud 63 Y. Therefore, we will pick the first value to avoid violating the slow roll 
condition. 

By adopting Planck, 2015 results [37], [62]-[63], we have ^ = 0.05Mpc“^, = 0.9682 ± 

0.0062, As = 2.196 x 10“^, and < 3.6 x 10“^Mpi, also by using ~ 100, Eqs.(4.12)- 
(4.13) become 


1 3Wj;q\i 


61.1 — Nk — In 


Hk 


(4.14) 


Treh = exp 


^ (1 + 'Wreh ) ■V'eh 






(4.15) 


The equation of state parameter at reheating era, iCreh? for general inflationary potential 
is usually taken in the interval, —1/3 < iCreh < 1 [13]. But for Starobinsky inflation {R^- 
inflation), the interval that fits well with Planck, 2015 results is 0 < rcreh <1/3 [62]. 


4.2 Constraining reheating parameters by the present PMF in i?^-inflation 

The next step is to relate the reheating parameters to the present PMF which is constrained 
by the results of Planck, 2015 [9]. As the conformal invariance of electromagnetic field is 
restored after inflation, the present super Hubble magnetic field Bq is redshifted since the 
end of inflation [59] as. 




^end 

(1 + ^end)^ 


(4.16) 


where, Hgnd is the magnetic field at the end of inflation and Zend is the redshift at the end of 
inflation. Hence, at the end of inflation we can write. 


CiQ CiQ tteq CLj-qIi 

1 + = - =- 

^end ^eq ^reh ^end 


-f^rad \ Pj J 


(4.17) 


where, (5.7 X is energy density of radiation today. 

We substitute (4.17) into (4.16) and make sure there is no backreaction problem at the 
end of inflation, ^ Pend- Since, ps — BA then combining (4.16) and (4.17) yields the 
constraint in reheating parameters from the present PMF [59], 


^ 1/2 

The upper limit of present PMF calculated by Planck, 2015 is ^ 10“^G. Therefore, the lower 
limit of reheating parameter is, i?rad > 1-761 x 10“^. But this limit is independent of the 
inflationary model. Hence, for a more accurate limit of i?rad associated with inflation, 
one can use the result of [49] and adopting the Planck, 2015 results, at which the middle 
value, ^ 58.5, if rcreh is not constant [37]. It implies that, i^^ad > 6.888, which is three 
orders of magnitude more than the previous one. Therefore, the reheating in this case is 
more constraint by the inflationary model bound. 

In order to constrain the reheating energy scale, one can substitute (4.18) into (4.4) to 
obtain. 
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l-3^reh 


6(1+Wreh) 


Preh ^ Pend 



(4.19) 


Therefore, one has to find pend foi* fhe model of inflation. Also, the lower limit of reheating 


energy density should be greater than BBN energy, > Pnuc, where pnuc is in the order of 
lOMeV). 

We can use the upper bound of the energy density of inflation derived from WMAP7, 
(Pend)cMB ^ 2.789 X [56] and the lower limit of i?rad(> 6 .888) and substitute them 

into (4.17), one can find the upper bound of, 1 + Zend < 6.828 x 10^^. If we substitute 
the redshift into 4.16, we can find the upper limit of magnetic field at the end of inflation, 
Send < 4.662 X 10^^G(2.3541 x 10“^Mpj), where we have used, IG 3.3 x 10“^^Mpp Hence, 
the energy density of the magnetic field is {pBend)cMB ^ 1-184 x 10“^^Mpp It is ten or¬ 
ders of magnitude more than the upper limit of pend found from CMB [56]. Therefore, the 
backreaction problem can be avoided easily. 

On the other hand, we can use Eq.(4.11) to find the energy density at the end of 
inflation (eiy l,Pend = fiend)- By adopting the value of A(^ 4.0 x 10“^Mpi), calculated 
from the amplitude of the CMB anisotropies [30], we find (Pend)R2_inflation — B1 x 10“^^Mpp 
Similarly, substituting of this value with the limit, i?rad(> 6.888) into (4.16)-(4.17) implies 
that, 1 + 2;end < 5.414 x 10^^ and (P5end)R2 -inflation ^ 4.6788 x 10“^^Mpp which is free more 
from backreaction problem. Therefore, (PBend)R2-inflation calculated by inflationary model 
and the present limits of PMF found by Planck, 2015 puts more constraints on the PBend 
than the constraints found by (pSend)cMB with present limits of PMF. 

Similarly, one can plot both A^^eh ^^nd T^eh as a function of rig by using Eq.(4.15), 
see Fig. 11-12. As we can see from Fig. 12, all curves of possible intersect at Treh ~ 
4.32 X lO^^GeV and ris ^ 0.9674, the value of spectral index fits well with the range of 
Planck, 2015 for i?^-inflation. This temperature is much more than the range of reheating 
temperature obtained from CMB in the context of LFI, see Eq.(4.11) of Ref. [56]. As we use 
the upper limit of in Eq.(4.14), we can consider the above value as the upper limit of 


Treh' 


Thus, adopting the this temperature for all Wj-qy^ models of reheating, enables us to 
constrain the A^reh on the range 1 < A^reh < 8.3, for all possible values of ,idi.eh see Fig. 13. 
The average value of A^^eh is relatively low, A^^eh — 4.7. Also, from, (4.6) one can find, 
Preh ~ 3.259 X 10“^^Mpj, which is in the same order of magnitude of the upper limit of 
(Pend)R2_inflation ^wo Orders of magnitude less than the upper limit of (pBend)cMB* Iciter- 
estingly enough this value of pjg^(1.035 x lO^^GeV) is one order of magnitude less than the 
range of Pend? obtained from CMB for the large field inflation [56], see Eq.(4.1) of the same 
reference. Therefore, the instantaneous reheating at which, pend ~ Preh? can be manifested 
by this result. 

At the end of this section, we constrain the upper limit of reheating energy density, 
Preh? based on the lower limits of both A^reh and reheating parameter, i?rad- The last one 
was obtained from the present upper limit of PMF by Planck, 2015. Hence, Eq.(4.1) can be 
written as. 



(4.20) 
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Figure 11. The plot of reheating e-folds, A^reh, versus spectral index, Ug, at the end of infiat ion, 
for some values in —0.3 < rCreh < 1- They all intersect at ~ 0.9674, which lies well in the Planck, 
2015 range, = 0.9682 ± 0.0062 


Treh (GeV) 



Figure 12. The plot of reheating temperature, Treh, versus spectral index, Ug, at the end of R^- 
infiation, for some values in —0.3 < Wj-eh < 1- They all intersect into Treh ^ 4.32 x 10^^GeV at 
Ug ~ 0.9674. The value of the temperature is below (pend)^^^ and the spectral index lie in the Planck, 
2015 range, = 0.9682 ± 0.0062. 


If we use (pend)cMB’ Upper limit, (preh)cMB ^ 8.480 X 10 which is much more than 

the lower limit derived from WMAP7 for both large and small field inflation [56]. However, 
if (/?end)R2 .inflation i^ USed, the Upper limit is (preh)R2 .inflation < 3-344 X It is still 

much more than the lower limit derived from WMAP7 for both large and small held inhation. 
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Figure 13. The number of e-folds, A^reh, during reheating, for all possible values of the equation 
of state parameter, —1/3 < Wreh < 1, at rig ~ 0.9637. The range of e-folds is, 1 < A^reh < 8.3. 
In inflat ion, the values of state parameter lies onto 0 < Wj-eh < 1/3, hence the middle value is 
A^reh ^ 4.7. 


5 Constraining the present PMF from the magnetogensis 


In this section, we will constrain the value of present PMF, based on both the predicted scale 
invariant magnetic field generated in inflationary era by f^FFmodel in i?^-inflation, which is 
discussed in section. 3, and by using the constraint values of reheating parameters, which are 
calculated in section. 4. In order to do that, we need to impose some necessary assumptions. 
First as shown in [24], the magnetic field energy density, pn, scales as 1/a^ independently 
of the dominant constituent in the reheating era and the eras come after. That is basically 
the implication of Eq.(4.16). Second, and to insure that, this model does not suffer from the 
backreaction problem, we generalize the validity of Eq.(102) in [24] to the end of inflation. 


dpE dpB 

dhik dink 


^ Pend') 

end 


(5.1) 


where, the EM spectra in (5.1) are calculated by (3.10)-(3.11). The third assumption is that, 
the reheating is not going to affect the magnitude and the shape of EM spectra. 

As a result of the above assumptions, one can neglect the electric field in Fig. 9 without 
specifying the constraints of electric conductivity of the reheating (ctc ^ H) which may lead 
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to zero electric field and constant magnetic field. Also, because of relatively small value of 


Areh? we can assume that, 


dpB 


dink 


dpB 


end 


d\nk 


reh 


Since, we use Planck units in computation 


(Mpi = 1), then the scale of the spectra is in (Mp^). By using the upper limit of Hubble 
parameter during inf 
from Fig. 9-10 is 


ation, H — 3.6 x 10 ^Mpi, the magnitude of magnetic spectra obtained 
- 3.842 X 10 ^^Mpp This value is well below both the upper limit of 

end 

(/Oend)cMB (/»end)R2_i„flatioii- However, it is One order of magnitude more than ipB^^a)cMB 


d\nk 


and much more than 


'end / R2 —inflation * 


Also, 


dpB 


d\nk 


end 


^ (Preh)R2_inflation’(/^reh)cMB* 


last result shows that the inflationary magnetogensis may play significant role during the 
reheating era. 

By taking the magnitude of the magnetic spectra as equivalent to fhen the upper 

bound of magnetic field at the end of inflation, Send < 3.76 x lO^^G. Therefore, if we use the 
upper bound of redshift derived from (pend)cMB’ ^ + ^end < 6.828 x 10^^, the upper limit of 
the present PMF is, Bq < 8.058 x 10“^G. It is in the same order of magnitude of the upper 
bound of PMF obtained by Planck, 2015. However, if we use the upper bound of redshift 
derived from (Pend)R2_inflation’ 1 + ^end < 5.414 X 10^^, the upper limit of the present PMF 
will be Bq < 1.282 x 10“^G. That is even higher than the galactic or interplanetary magnetic 
field which is of the order of 10“^G. Therefore, the second limit is too weak. 


6 Summary 

In this paper we have shown that the scale invariant PMF can be generated by the simple 
model pFF^ Eq.(l.l), in i?^-inflationary model, which is mostly favored by the latest result 
of Planck, 2015 [37]. Unlike, in non-standard inflationary models like NI [39], LFI [40], and 
some standard models of inflation [24] at which pFF model suffers from the backreaction 
problem, we can avoid this problem in i?^-inflation. It is easily to avoid this problem as long 
as, the rate of inflationary expansion, i7, is in the order of or less than the upper bound 
reported by Planck (< 3.6 x 10“^Mpi) [37]. In principle, the electric spectrum starts exceeds 
the magnetic spectrum at relatively high, H{> 0.2Mpi). The corresponding e-folds number 
is, N ^ 68. Thus, we can consider these two values as upper bounds to the model. We 
do this investigation for both simple exponential (de Sitter) and power law expansion. At 
sufficiently high e-folds number, N^ there is no significant differences in their results. 

The second main result is constraining the reheating parameters from the upper limits 
of PMF reported by Planck, 2015 [9]. In the context of i?^-inflation, we calculate the lower 
limits of the reheating parameter, i?rad > 6.888. Also, we find the other reheating parameters 
based on the upper limit of energy density at the end of inflation calculated from CBM 
data, (/Oend)cMB < 2-789 x and from the inflationary model, (Peiid)R2_inflation - 

1.1 X 10“^®Mpj. As a result, we find that the magnetic field energy density at the end of 
inflation as (ps,„d)cMB < 1-184 x and (pSend)R2-inflation < 4-6788 x^lO-^^M^p 

Similarly, we find the upper limit at the end of reheating, (preh)cMB ^ 8.480 x 10“^^Mpj and 
(Preh)R2_ inflat ion ^ 3.344 X 10“^^Mpp All of foregoing results are more than the lower limit 
derived from WMAP7 for both large and small field inflation [56]. These results show the 
significance of PMF during reheating era. 

On the other hand, we constrain the reheating parameters by using the Planck infla¬ 
tionary constraints, 2015 [37] in the context of i?^-inflation. The upper limit of reheating 
temperature and energy density for all possible values of ,R;reh are respectively constrained 
as, Treh < 4.32 x lO^^GeV and preh < 3.259 x 10“^^Mpj at 0.9674. This value of 
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spectral index is well consistent with Planck, 2015 results. Adopting T^-eh for all models 
of reheating, enables us to constrain the A^^eh on the range 1 < A^eh < 8.3, for all possible 
values of ,r(;reh- 

At the end, we constrain Sq, from the scale invariant PMF, generated by f^FF in 
i?^-inflation in section. 3. From the PMF spectra. Fig. 9-2.8, we find that the upper limit of 
magnetic field in the end of inflation is Send < 3.76 x lO^^G. Therefore, if we use the upper 
bound of redshift derived from (pend)cMB’ fhon Bq < 8.058 x 10“^G. It is in the same order 
of PMF obtained by Planck, 2015. However, if we use the upper bound of redshift derived 
from (Pend)R2_inflation’ f^on Bq < 1.282 X 10“^G. That is even higher than the interplanetary 
or galactic magnetic field which is of the order of 10“^G. Therefore, the second limit is too 
weak. 

In order to achieve the scale invariant PMF by this model, the free index of the coupling 
function has a relatively high values, a ^ {—7.44, 7.44}. However, at a = 2, which is the 
typical value, we cannot generate scale invariant magnetic field although we can avoid the 
backreaction problem. The main problem with this model is the value, \a\ ~ 7.44, which 
is out of the dynamo constraints imposed by CMB, BBN, and Faraday rotation, RM, see 
Eq.(94) of [24], see Fig.3 of the same reference. In fact those limits are derived mainly for 
exponential and large field inflation models. Therefore, this subject needs more investigation 
on the context of i?^-inflation. 

If the above problem, the relatively high value of a, is not resolved and i?^-inflation 
shows more consistency with cosmological observations, the results of this research would 
go against the simple inflation model f^FF^ as way of generating PMF. These difficulties 
add constraints to this model, in addition to those found by other researches, such as new 
stringent upper limits on the PMF, derived from analyzing the expected imprint of PMF on 
the CMB power spectra [66], bi-spectra [67], tri-spectra [68], anisotropies and B-modes [69], 
and the curvature perturbation and scale of inflation [69]-[70]. However, if there is a way to 
justify the relatively high a and the approximation of quasi-const ant A, we assume to derive 
Eq.3.16, then pFF model may be viable in the context of i?^-inflation. Nevertheless, the 
agreement between the result of this model and the upper limit of present PMF found by 
Planck, 2015 cast some credits to this model. 

One way to extend this research, is to consider calculating 3-point statistics such as 
{BBB)^ {BB(j)) or {BBC) [71]-[73]. That may enable us to investigate the footprint of the 
PMF on CMB. It might be addressed in a future research. 
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